Widely tunable single photon source 
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Abstract: We theoretically and experimentally investigate the spec- 
tral tunability and purity of photon pairs generated from spontaneous 
parametric down conversion in periodically poled KT1OPO4 crystal with 
group- velocity matching condition. The numerical simulation predicts that 
the purity of joint spectral intensity (Pjsi) and the purity of joint spectral 
amplitude (Pjsa) can be kept higher than 0.98 and 0.81, respectively, when 
the wavelength is tuned from 1460 nm to 1675 nm, which covers the S-, 
C-, L-, and U-band in telecommunication wavelengths. We also directly 
measured the joint spectral intensity at 1565 nm, 1584 nm and 1565 nm, 
yielding Pjsi of 0.989, 0.983 and 0.958, respectively. Such a photon source 
is useful for quantum information and communication systems. 
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1. Introduction 

As it is well known, classical widely tunable laser sources play an important role in many fields, 
for example, spectroscopy, remote sensing, metrology, and optical communications (HEI. How- 
ever, the quantum counterpart of the tunable laser, the widely tunable single photon source, 
have been less investigated. Wavelength tunable single photon sources have been generated 
from quantum dots, but the tunable ranges achieved until now have been narrow EH). 

In this paper, we present a widely tunable single photon source based on spontaneous para- 
metric down conversion (SPDC) in a periodically poled KT1OPO4 (PPKTP) crystal with group- 
velocity matching condition. The concept of group- velocity matching in SPDC was introduced 
by Keller and Rubin 0, and by Grice and Walmsley O in 1997. PPKTP crystal with group- 
velocity matching condition was experimentally demonstrated by Konig and Wong Q for 
second-harmonic generation in 2004. Later, this condition was applied in SPDC for intrinsi- 
cally pure photon state generation. Under this condition, the signal and idler photons from 
SPDC have no spectral correlation, i.e., they are intrinsically pure. Therefore, there is no need 
to employ bandpass filters to obtain high-purity heralded single photons. Consequently, such 
photon sources are much brighter than the traditional sources and have showed high brightness 
in experiments, as reported by Evans, et al, for entangled photons H); by Eckstein, et al, for 
two-mode squeezer (9); and by Yabuno, et al, for four-photon interferometry [flOl . 

Besides the high purity and high brightness, in this research, we focus on another important 
merit of this source: the wide spectral tunability. From the numerical calculation, we found 
that the purity of joint spectral intensity (Pjsi) and the purity of joint spectral amplitude (Pjsa) 
can be kept higher than 0.98 and 0.81, respectively, as the wavelength is tuned from 1460 nm 
to 1675 nm. To verify this simulation, we directly measured the joint spectral intensity, and 
obtained high P JSI of 0.989, 0.983 and 0.958 at 1565 nm, 1584 nm and 1565 nm, respectively. 

2. Theory 

In the process of SPDC, one pump photon is split into two lower energy photons, the signal and 
the idler. The two-photon component in SPDC can be expressed as ifTTIl 



where /(fi) s ,cty) = (j)(cQ s ,CQi)a(cQ s + ooi) is the joint spectral amplitude (JSA), (j)(co S: (Oi) and 
a((D s + COi) are the phase matching amplitude and the pump envelope amplitude. CO is the 
angular frequency, is the creation operator and the subscripts s and i denote the signal 
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Fig. 1. Examples of (a) pump envelope intensity (PEI), (b) phase matching intensity (PMI) 
and (c) joint spectral intensity (JSI). (c) is the product of (a) and (b). The width of PEI 
is proportional to the bandwidth of the pump, while the tilting angle of PEI is fixed, 45 
degree to the minus direction of the horizontal axis. The width of PMI is proportional to 
the inverse of the crystal length. The direction of PMI is not fixed, but determined by the 
group velocities of the signal, the idler and the pump. When the width of PEI and PMI are 
equal and the PMI is at 45 degree to the positive direction of the horizontal axis, JSI can 
achieve a subcircular shape and maximum purity. 



and the idler photons, respectively. Assuming the spectrum of the pump laser has a Gaus- 
sian distribution with a bandwidth of o p , the pump envelope intensity can be written as 
\cc(co s + C0i)\ 2 = exp[— ( ^— ^ COp ) 2 ]. Under the collinear condition, the phase matching inten- 



sity can be written in the form of \(j)(co S: (0i)\ 2 = [sinc(^)] 2 , where Ak = k p — k s — k[ — ^ 
is the difference between the wave vector of the pump (k p ), the signal (k s ), the idler (kf) and 
an extra vector (^) introduced by periodical poling of the crystal. L and A are the length and 
poling period of the SPDC crystal. Figure fl]a-c) shows examples of the pump envelope in- 
tensity \a(co s + C0i)\ 2 , phase matching intensity \(j)((Q s , (Qi)\ 2 and joint spectral intensity (JSI) 

\f(C0 s ,C0i)\ 2 . 

Considering the Taylor expansion of the three wave vectors k^ = ^(co^o) +^(®juo)(^ju — 
tftyio) -\-o{{(Dn — CQ[io) 2 )> (M = PjS,i) 9 the difference of them can be approximated as Ak = 
A&(°) + Ak^ + .... In our experiment, the phase matching condition is satisfied at co s = C0i = 
(Do = ^ and thus the zero order contribution disappears: 

2tl 

A*(°) = k p (2coo) - k s (coo) - ki(coo) - — = 0. (2) 

The first order term can be written as 

Ak^ = k f p (2coo)(co p - 2(Oo) - K(coo)(co s -(Dq)- ^(coo)(^ - <flo) (3) 
= Wpilcoo) -*£(fflb)](fi>5 - coo) + Wp^coo) -*-(fflb)](fi^ - «o), 

where "/" represents the first order differentiation with respect to the angular frequency. 

As analyzed in fT2lh the slope of the phase matching condition is determined by the gradient 
of Ak. 

gradA£= (f^,f^) = (V^fa)-V£fa),V^ (4) 

where the group velocity V g ^ is defined as V g ^ = = p ^ , (/i = p,s,i). From the view- 

point of the tilting angle, 0, 



where is the angle between positive direction of horizontal axis and the ridge direction of the 
phase matching intensity, as shown in Fig.[TJb). 

The purity, a parameter which describes the spectral correlation of the photon source, is 
defined as p = Tr(pj) = Tr(pf), where p s = Tri{\\\f S i) {ty s i\) is the reduced density operator 
of the signal, and Tr ; represents the partial trace on the subsystem i. This purity is determined 
by the factorability of the JSA, f(co Sj COi), and can be calculated numerically using Schmidt 
decomposition (THE). By applying Schmidt decomposition on the JSA, we can obtain 

f((D U Q> 1 )=T, j Cj<l)j(CQ 1 )(pj(Q>2), (6) 

where (t>j((D\) and (pj(a>i) are two orthogonal basis sets of spectral functions, known as the 
Schmidt modes, cj is a set of real, non-negative, normalized weighting coefficients with Ljcj = 

1. The Schmidt number is defined as K = — [ 14]. K indicates the number of Schmidt modes 

Lj c j 

existing in the two-photon state and thus can be viewed as an indicator of entanglement. The 
purity of the signal or the idler state equals to the inverse of the Schmidt number K: 

P = PjSA = Y, c4 j' ( 7 ) 

j 

In experiment, it is difficult to directly measure f(cO\ , CO2), the JSA. What we usually measured 
is \f(d) s , C0i) | 2 , the JSI, as have been widely demonstrated in previous experiments ll8ll9l fT0l[TT1l . 
To analyze such experimental data in experiment, the Schmidt decomposition can be applied 
on |/0„ CO;) | 2 , 

\f(C0 U CO2)\ 2 = Y,j^j(^l)(Pj(0>2)- (8) 

We define the purity of the joint spectral intensity as 

>'/.s/-I>/. (9) 

Pjsi is useful in experiment to characterize the joint spectral intensity. In this paper we focus on 
the parameter Pj$i. 

3. Numerical simulation 

With the Sellmeier equation in [7 1, we calculated Pjsi and Pjsa versus the wavelength, as shown 
in Fig. [2 The horizontal axis is the wavelength of the signal (or the idler, in degenerate case), 
which is two times of the pump wavelength. The bandwidth of the pump is set to match the 
length of the crystal, so as to achieve the maximum purity at 1584 nm. The optimal bandwidths 
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Fig. 2. Numerical simulation of the Pjsi and the Pjsa as functions of the wavelength. 



(a) Pjsi=0.987, at 1500 nm (b) Pjsi=0.992, at 1550 nm (c) Pjsi=0.992, at 1600 nm (d) Pjsi=0.988, at 1650 nm 
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Fig. 3. Numerical simulation of the joint spectral intensity (JSI) at 1500 nm, 1550 nm, 
1600 nm and 1650 nm. 



for Pjsi and Pjsa are slightly different. It is noteworthy that the Pjsi and Pjsa can be kept 
between 0.983 and 0.993, and 0.815 and 0.821, respectively, from 1460 nm to 1675 nm. The 
Pjsa is smaller than the Pjsi, mainly caused by the side lobes in the phase matching amplitude 
(f) (co s , C0i) . The Pjsi and Pjsa are relatively low in the regions around 1450 nm and 1700 nm. This 
is because the same bandwidth of the pump is used for all the wavelengths, but this bandwidth 
is optimum only for 1584 nm. If the bandwidth of the pump is optimized for 1450 nm or 1700 
nm, we can also achieve a higher values for those regions. Figure[3a-d) shows several examples 
of the JSI at 1500 nm, 1550 nm, 1600 nm, 1650 nm, respectively. 

4. Experiment and results 

The experimental setup is shown in Fig. [4] Picosecond laser pulses (temporal duration ~ 2 ps, 
central wavelength was tunable from 700 nm to 1000 nm) from a mode-locked Titanium sap- 
phire laser (Coherent, Mira900) were used as the pump source for the SPDC. Pump pulses with 
power of 50 mW passed through a 30-mm-long PPKTP crystal with a poling period of 46. 1 fim 
for type-II SPDC. The temperature was maintained at 32 °C for the PPKTP crystal. The down- 
converted photons, i.e., the signal (vertically polarized) and the idler (horizontally polarized) 
were separated by a polarizing beam splitter (PBS) and collected into two single-mode fibers 
(SMF). Then the photons were filtered by two bandpass filters (BPFs, Optoquest), which had a 
filter function of Gaussian shape with an FWHM of 0.56 nm and a tunable central wavelength 
from 1560 nm to 1620 nm. Finally, all the collected photons were sent to two InGaAs avalanche 
photodiode (APD) detectors (ID210, idQuantique) connected to a coincidence counter (Ortec 
9353). 

To measure the JSI of the photon pairs, we scanned the central wavelength of the BPF1 and 
BPF2, and recorded the coincidence counts. BPF1 and BPF2 were moved in 0.1 nm per step 
and 60 by 60 steps in all. The coincidence counts were accumulated for 10 seconds for each 
point. With the pump wavelength set at 782.5 nm, 792 nm and 807.5 nm, we measured the JSI 
at 1565 nm, 1584 nm and 1615 nm, respectively. 

The measured Pjsi were 0.989, 0.983 and 0.958, respectively, as shown in Fig. [5] The max- 
imum coincidence counts rate decreased from 1075 cps at 1565 nm to 281 cps at 1615 nm, 
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Fig. 4. The experimental setup. LPF (long-wave pass filter), PBS (polarizing beam splitter), 
SMF (single mode fiber), BPF (bandpass filter), APD (avalanche photodiodes). 




(a) Pjsi=0.989, at 1565 nm (b) Pjsi=0.983, at 1584 nm (c) 
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Fig. 5. Experimentally measured joint spectral intensity (JSI) at 1565 nm, 1584 nm and 
1615 nm. 
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Table 1. Comparison of the theoretical, convolved and experimental Pjsi. 



due to the decrease of the quantum efficiency of the APDs. As depicted in Fig. [2 the theoret- 
ical Pjsi for 1565 nm, 1584 nm and 1615 nm, were 0.992, 0.992 and 0.991, respectively. In 
the simulation of the theoretical Pjsi, we assumed that the bandwidth of two BPFs was narrow 
enough. However, in experiment the bandwidth of our BPFs (0.56 nm) was comparable to the 
bandwidth of the signal and idler (1.1 nm). So, it was necessary to consider the convolution 
effect of the two BPFs. We repeated the simulation with the real filter function and obtained 
the convolved Pjsi as 0.995, 0.995 and 0.995. The three measured Pjsi were very close to these 
convolved values, as shown in Tab. Q] The small disparities may have been caused by the the 
dark counts in the APDs, or small differences in the BPFs. 

It was noteworthy to compare the theoretically expected bandwidth and experimentally meas- 
ured one. The measured spectrum was a convolution of the original spectrum and the filter 
function of the BPFs. The BPFs had a filter function of Gaussian shape with an FWHM of 
0.56 nm, while the spectra of the signal and idler were also in a Gaussian shape. Thus we can 

calculate the their convolution as FWHM con = yj 0.56 2 + FWHM t 2 he , where FWHM the was the 
theoretical calculated FWHM of the original JSI. The experimentally measured FWHM exp can 
be obtained from the marginal distribution of the data in Fig. [5j We compared the theoreti- 
cal bandwidth and the experimentally measured bandwidth in Tab. [2 The convolved FWHMs 
were consistent with the measured ones. The small difference may arise from that the FWHM 
of pump laser was slightly different from the optimal FWHM in the simulation. In addition, the 
FWHM of the pump laser and the FWHM of the two BPFs were slightly wavelength-dependent. 



5. Discussion and future 

We need to clarify that purity of joint spectral intensity, Pjsi, is different from the traditional 
definition of spectral purity p, or Pjsa> There is no information on phase of the two-photon wave 
packet in JSI. Any phase would reduce the purity of the heralded single photons by introducing 
temporal correlations. Thus, we cannot directly obtain purity of the heralded single photons 
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Table 2. Comparison of the convolved bandwidth and the experimentally measured band- 
width. 



from JSI. As shown in Fig. O Pjsa is usually smaller than Pjsi, due to the side lobes in the 
JSA. However, the detrimental effect of the side lobes is not severe, according to the previous 
experiments. Several groups have reported high- visibility interference between the signal and 
idler photons from such PPKTP crystals, e.g., visibility has achieved 94% in [8] and 98% in 
|[T0l . Especially, the interference visibility was 95% for both PPKTP crystal and custom-poled 
crystal in 11751 . 

To reduce the effect of side lobes in JSA, a new method of engineering the nonlinearity profile 
of the poled crystals has been proposed recently [fl"5l . With this method, the phase matching 
amplitude can be tailored to approximate a Gaussian function. Therefore, the side lobes can be 
significantly suppressed and the Pjsa can be improved from 0.81 to 0.99 at 1576 nm [TBI . They 
used higher-order poling to modulate the crystal. In this process, the material property of KTP 
was not changed. Therefore, in principle, widely tunable photon sources with higher purity 
could be generated with such custom-poled crystal. However, this method had the drawbacks 
of smaller effective nonlinearity and thus relatively lower photon pair rate [1131 compared to 
the traditional PPKTP crystals. Alternatively, the purity can be improved by using a wide-band 
BPF to filter out the side lobes, but this method slightly decreases the count rate. 

This spectrally pure photon source can be easily transformed into an entangled photon source 
by adding two calcite prisms as displacers 1U [16) or setting it into a Sagnac loop configuration 
ifTTl IT81 . This photon source can be used for pulsed two-mode squeezed state generation by 
combining the signal and idler in the degenerate condition f[T9l . Another possible application 
for this photon source is wavelength-multiplexing based multiparty quantum communication 
system. If this PPKTP crystal is pumped by a laser frequency comb, e.g., with wavelengths 
distributed from 730 nm to 835 nm, we can generate serial pure photons distributed from 1460 
nm to 1675 nm. These photons can be used for a wavelength-multiplexed system, which is one 
step for practical multiparty quantum communication systems. 

The classical broadly tunable laser source have showed a tremendous impact in many and 
diverse fields of science and technology HI . As the quantum counterpart of classical broadly 
tunable laser source, our widely tunable single photon source also has the potential to play 
an important role in quantum information and communication systems, when the wavelength 
tunability are necessary. 

6. Summary 

In summary, we have demonstrated the generation of widely tunable and highly pure photons 
from PPKTP crystal. In the theoretical simulation, we found the wavelength can be tuned from 
1460 nm to 1675 nm, with Pjsi over 0.98 and Pjsa over 0.81. In experiment, we measured 



the JSI at 1565 nm, 1584 nm and 1565 nm, and achieved P JSI of 0.989, 0.983 and 0.958, 
respectively. This result was well consistent with our theoretical simulation. We also discussed 
the future application of this source. 



Acknowledgments 

The authors are grateful to M. Yabuno, P. G. Evans and T. Gerrits for helpful discussions. This 
work was supported by the Founding Program for World-Leading Innovative R&D on Science 
and Technology (FIRST). 



